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ABSTRACT 

In this paper we present a study of chemical abundances in six star-forming regions. Stellar parameters and metallicities are 
derived using high-resolution, high S/N spectra of weak-line T-Tauri stars in each region. The results show that nearby star- 
forming regions have a very small abundance dispersion (only 0.033 dex in [Fe/If]). The average metallicity found is slightly 
below that of the Sun, although compatible with solar once the errors are taken into account. The derived abundances for Si and 
Ni show that the observed stars have the abundances typical of Galactic thin disk stars of the same metallicity. The impact of 
, these observations is briefly discussed in the context of the Galactic chemical evolution, local inter-stellar medium abundances, 
■ and in the origin of metal-rich stars in the solar neighbourhood (namely, stars more likely to harbour planets). The implication 
' for future planet-search programmes around very young, nearby stars is also discussed. 

' Key words, planetary systems: formation - ISM: abundances Stars; abundances - Galaxy: abundances - solar neighbourhood - 
, Stars: formation 



1. Introduction 

As of November 2007, about 250 extra-solar planets 
have been discovered orbiting solar-type star^j, most 
of which were detected thanks to the development of 
hig h precision radial-vel ocity instruments (for a review 
see Udry fc Santos! [20071 ). Complementary high accuracy 



spectroscopic studies have shown that stars hosting giant 
planets are particul arly metal-rich when compared with 
single field stars (e .g. Gonzale j ll997t Gonzalez et a.1 ]l200ll 



Santos et all I200lll2004 , 2003 : iFischer fc Valentil l2005f) . 



* Based on observations collected at Paranal Observatory, 
ESO (Chile) with the UVES and flames/uves spectrographs 
at the VLT-Kueyen spectrograph (run IDs 075.C-0272(A) and 
076.C-0524(A), respectively), with the feros spectrograph at 
the eso/mpi 2.2-m telescope (program ID 070.C-0507(A)), as 
well as with the SOPHIE spectrograph at the OHP observatory, 
France. 

^ For a continuously updated list see tables at 



This observation is helping to better understand the pro- 
cesses of planet formation and evolution. 

The high metal content of some of the planet-hosts 
has raised a number of interesting questions regarding 
their origin. It has been proposed that the high metallici- 
ties could be the result of the in-fall of planetary (metal- 
rich) material into the stellar outer conv ective envelope 
(e.g. .Gonzalez 1998t Murray et al. 19981 ). although cur 



rent r esults do not support this hyp othesis (jSantos et al 



2004L l2005t iFischer fc Valentil l2005h . The excess metal 
licity observed likely reflects the higher (average) metal 
abundances of the clouds of gas and dust that originated 
the star and planet systems. Such conclusions are also 
supported by the most recent pl anet formation models 



(|lda fc Linll2004 iBenz et al.ll2006l ) 



Because the large majority of planet hosting 
stars are solar-neig hbourhood thin-disk objects (e.g. 



Ecuvillon et al.l 120071 ) . one would expect nearby T-Tauri 



http://www.exoplanets.eul or http://www.exoplanet.eu 



stars and star-forming regions [sfrs] that are metal-rich 
to exist. Their existence is expected if a well-defined 
age-metallicity relation exists in the solar-neighbourhood 
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2004 iHavwoodl l2006l) . Previous studies 
James et al.l2006( ). however, have not suc- 



ceeded in finding evidence of high metalhcity, T-Tauri 
stars in nearby SFRs. Whether this is due to the absence 
of such regions, or to the relatively small number of stars 
and SFRS observed was not clear. 

Considering the strong metallicity-giant planet con- 
nection, the detection of nearby, metal-rich SFRs would 
provide preferential targets for future planet searches 
around young T-Tauri stars. The detection of planets 
orbiting T-Tauri stars would provide observational con- 
straints of paramount importance for the study of planet 
formation scenarios, since it gives us some hint about the 
timescales for planetary formation. 

In thi s paper, we continu e the metallicity programme 
started in James et al.l ( 20061 ). by increasing the number of 
observed young T-Tauri stars in each SFR, and by adding 
abundance determinations for objects that have not yet 
been studied. In all, we present accurate stellar parameters 
and chemical abundances for 19 weak line T-Tauri stars 
[WTTSs] in six nearby SFRs, namely Chamaeleon, Corona 
Australis, Lupus, Rho-Ophiuchus, Taurus and the Orion 
Nebula Cloud. The results of this study are presented and 
their implications discussed. 

2. The sample 

For WTTSs in the Chamaeleon, Corona Australis, Lupus, 
Rho-Ophiuchus, and Taurus SFRs, the sample used in the 
current paper to derive stellar parameters and chemical 
abundances was chosen in an identical manner to that in 



James et al.l (|2006f ) , to which we refer the reader for exten- 
sive details. In brief, wtts catalogues for each sfr were 
constructed from ROSAT All-Sky Survey [rass] detections 
in and around each sfr. We then selected those stars hav- 
ing spectral types between G and early-K and visual mag- 
nitudes of 12 or brighter, which also have weak signatures 
of infrared excess and magnetic activity {e.g., having Ha 
equivalent widths <<10A). We also restricted our sample 
to objects exhibiting substantial equivalent widths [ews] 
of the lithium resonance line at 6708A (whose substantial 
abundance is indicative of youth), as well as those stars 
which arc not components of multiple systems, as judged 
from existing kinematic d ata. 

In I James et al. ( 20061) . we presented the results of the 
analysis of a sub-sample of these stars, whose spectra were 
obtained using the feros spectrograph (2.2-m eso/mpi 
telescope). The data were used to study the chemical 
abundances of wttss in 3 southern star-forming regions 
(Chamaeleon, CrA and Lupus). Unfortunately, the small 
aperture of the telescope did not allow us to obtain high 
S/N data for that many interesting targets. The number 
of stars observed was thus quite low, and the S/N of their 
spectra was often not of a quality sufficient to derive accu- 
rate stellar parameters and abundances. The use of a high 
resolution spectrograph installed on a larger telescope was 
thus mandatory to improve our results. In t his paper, we 
present an extension of the study initiated in I James et al 
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Fig. 1. Sample UVES spectrum for one of the stars 
in our sample (RX J1159.7— 7601) in the region around 
5950A. The spectrum presented has a S/N ratio of ~140. 
Gaussian fits to most of the lines in the spectral domain 
presented are shown (dashed line). The EWs of the lines 
go from 12 to 130 mA, the typical values for the Fei lines 
in WTTS spectra used to derive stellar parameters and 
metallicity. 



( 20061 ). using new high S/N, high resolution spectroscopic 
observations, predominantly obtained using 8m-class tele- 
scopes. 

After our new stellar spectra were reduced (see Sect. [3]) 
we restricted a detailed spectroscopic analysis to narrow- 
lined stars, i.e., those having low projected equatorial ro- 
tational velocity (<~ 15kms~^). Calculations of projected 
equatorial rotational velocity [v sin i] for all stars were de- 
rived using the width of the cross-correla t ion fu nction, us- 
ing a method des cribed by Melo et al] (|2001l ) (see also 
James et al.l 120061 ) . This is imperative because higher ro- 
tational velocities pose important limitations on the ac- 
curacy of derived stellar parameters and chemical abun- 
dances (see Sect.lD). 

One-dimensional kinematic data (radial velocities 
[rvs] ) of the target stars were compared with the expected 
systemic velocities of their parent sfr (using values cal- 
culated and compared in James et al. 2006). Target stars 
were only considered for a detailed spectroscopic abun- 
dance analysis if their rvs were invariant (i.e., not in a bi- 
nary system) and they matched that of their natal cloud. 
A full description of the membership, mass and age prop- 
erties of our target stars, as well as details on the rapid 
rotators, is presented in James et al. (in prep.). 

For stars in the Orion Nebula Cloud, data and 
membership information were taken from the studies 
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Table 1. List of WTTS observed in each SFR and derived stellar paranreters, nretallicities, and estimated v sini values. 
The number of Fe i and Fe ii lines used for the spectroscopic analysis is also listed, together with the RMS around the 
aver age abundance for e ach set of lines. Stellar parameters and [Fe/H] derived from feros spectra ([4]) were taken 



from I James et al.l (120061 ) 



Star 


Toff 


log gspec 




[Fc/H] 


Insfl" 


N 


a 


V sin i 




[K] 


[cms^'^] 


[kms^^] 






(Fe l,ll) 


(Fe I, ll) 


[kms^^] 


Cha 


















RX J1140.3-8321 


4664± 94 


4.17±0.52 


2.11±0.15 


-0.22±0.11 


[1] 


31,4 


0.11,0.29 


11 


RX J1140.3-8321 


4814±154 


4.07±0.92 


2.15±0.23 


-0.14±0.15 


[4] 








RX J1158.5-7754A 


4774±101 


4.08±0.41 


2.11±0.15 


-0.28±0.11 


[1] 


29,5 


0.10,0.21 


11 


RX J1158.5-7754A 


4810±142 


3.87±0.51 


2.11±0.22 


-0.26±0.15 


[4] 








RX J1159.7-7601 


4740± 81 


4.37±0.33 


1.70±0.12 


-0.08±0.08 


[1] 


32,5 


0.08,0.17 


9 


RX J1159.7-7601 


4836±146 


4.23±0.33 


2.07±0.22 


-0.18±0.16 


[4] 








RX J1233.5-7523 


5524± 53 


4.49±0.16 


1.22±0.07 


-0.03±0.06 


[1] 


35,11 


0.06,0.07 


6 


RX J1233.5-7523 


5469± 62 


4.40±0.22 


1.52±0.10 


-0.09±0.08 


[4] 








CrA 


















CrAPMS2 


5400±164 


4.06±0.45 


2.41±0.30 


-0.09±0.18 


[4] 








CrAPMS-4SE 


5236± 58 


4.33±0.17 


2.15±0.09 


-0.02±0.08 


[1] 


35,8 


0.07,0.08 


10 


CrAPMS-4SE 


5335±101 


4.59±0.41 


2.24±0.17 


-0.02±0.12 


[4] 








RX J1839.0-3726 


5150± 66 


4.40±0.22 


2.45±0.11 


-0.10±0.08 


[1] 


36,9 


0.07,0.11 


10 


RX J1839.0-3726 


5293±127 


4.51±0.42 


2.19±0.21 


-0.01±0.15 


[4] 








Lup 


















RX J1507.2-3505 


5169± 63 


4.31±0.33 


2.09±0.09 


-0.06±0.08 


[1] 


37,9 


0.07,0.17 


13 


RX J1507.2-3505 


5155±148 


4.37±0.49 


2.53±0.26 


-0.08±0.14 


[4] 








RX J1546.6-3618 


5223± 83 


4.51±0.22 


2.22±0.12 


-0.04±0.09 


[1] 


34,8 


0.09,0.10 


6 


RX J1546.6-3618 


5062±109 


4.09±0.53 


2.12±0.16 


-0.12±0.14 


[4] 








RX J1547.6-4018 


5006± 61 


4.29±0.24 


1.87±0.08 


-0.06±0.07 


[1] 


34,7 


0.07,0.12 


12 


RX J1547.6-4018 


5045±120 


4.18±0.61 


2.23±0.20 


-0.14±0.13 


[4] 








RX J1601. 1-3320 


5479±152 


3.70±0.35 


2.49±0.26 


-0.05±0.18 


[4] 








RX J1523.4-4055 


4923± 87 


4.74±0.37 


2.05±0.16 


-0.03±0.09 


[1] 


34,6 


0.09,0.19 


10 


Rho-Ohp 


















RX J1620. 7-2348 


4867±107 


4.40±0.29 


2.33±0.17 


-0.08±0.12 


[1] 


35,7 


0.12,0.13 


9 


ONC 


















JW157 


4877± 83 


3.57±0.28 


1.97±0.09 


-0.06±0.11 


[2] 


36,10 


0.10,0.14 


7 


JW589 


4927±108 


3.90±0.30 


1.95±0.13 


-0.12±0.12 


[2] 


33,5 


0.11,0.15 


14 


JW706 


4599± 79 


3.86±0.45 


1.50±0.10 


-0.17±0.09 


[2] 


30,4 


0.09,0.25 


12 


Tau 


















1RXSJ052146.7 


4921± 59 


4.05±0.29 


1.92±0.07 


-0.07±0.07 


[3] 


34,5 


0.07,0.15 


13 


1RXSJ053650.0 


4918±112 


4.47±0.29 


1.96±0.17 


-0.18±0.10 


[3] 


28,3 


0.09,0.12 


14 


1RXSJ053931.0 


5430± 53 


4.09±0.23 


1.87±0.08 


0.05±0.07 


[31 


36,12 


0.06,0.10 


11 



t Instruments used: [1] UVES: [2] FLAMES: [3] SOPHIE: [4] feros data from I James et al.1 (fioOd ) 



of ISicilia-Aguilar et al.l (|2005l ). iRhode et all (|200lf) . and 
Jones fc Walker] (jigSST ). 

We note that one star from the Chamaeleon SFR 
(RX J1201.7-7959 ) whose metallicity was derived in 
James et al.l (|2006f ) is now excluded from the sample, be- 
cause its Rvs indicate that it is a spectroscopic binary 
system. 

3. Observations 

Between April and September 2005, we obtained a se- 
ries of high resolution spectra for wttss in 4 young 
SFRs (Chamaeleon, Corona Australis, Lupus and Rho- 
Ophiuchus) using the uVES spectrograph on the vlt's 
Kueyen 8.2-m telescope (eso program ID 075.C-0272). 
The spectra were obtained in Service Mode using the 
Red 580 mode, and cover the wavelength domain between 



4800 and 7000A, with a gap between 5700 and 5850A. A 
binning of 2x2 pixels and a slit width of 0.9 arcsec was 
used, leading to a resolution of R=A/AA ^^50000. A sam- 
ple UVES spectrum used is presented in Fig.[T] 

To complement the uvES observations, six wttss 
in two more star-forming regions were observed. 
Observations of objects in the Orion Nebula Cloud [ONC] 
were carried out using the fibre link of the flames-UVES 
spectrograph (eso program ID 076.C-0524). The observa- 
tions were carried out in Service Mode between October 
2005 and March 2006. The CCD was used in a 1x1 bin- 
ning mode, and the resultant spectra have a resolution 
of ~50 000 and cover the same spectral domain as the 
UVES data. Finally, spectra for wtts in the northern 
Taurus SFR were obtained in December 2006 using the SO- 
PHIE fibre- fed spectrograph at the 1.93-m telescope of the 
Observatoire de Haute Provence [ohp], France. The spec- 
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Table 2. Abundances of Si and Ni for the stars in Tabled] that were observed with UVES, flames or SOPHIE spectro- 
graphs. The [Fe/H] values listed in Tabled] are also presented for comparison. 



Star 


[Fe/H] 


[Si/H] 


a([Si/H]) 


N(Si) 


[Ni/H] 


a([Ni/H]) 


N(Ni) 


Cha 
















RX J1140.3-8321 


-0.22 


-0.19 


0.08 


4 


-0.30 


0.13 


15 


RX J1158.5-7754A 


-0.28 


-0.27 


0.07 


4 


-0.34 


0.14 


14 


RX J1159.7-7601 


-0.08 


-0.05 


0.07 


5 


-0.13 


0.09 


18 


RX J1233.5-7523 


-0.03 


0.00 


0.04 


9 


-0.07 


0.04 


18 


CrA 
















CrAPMS-4SE 


-0.02 


-0.01 


0.07 


9 


-0.08 


0.13 


18 


RX J1839.0-3726 


-0.10 


-0.10 


0.10 


8 


-0.22 


0.11 


16 


Lup 
















RX J1507.2-3505 


-0.06 


-0.02 


0.09 


8 


-0.18 


0.08 


14 


RX J1523.4-4055 


-0.03 


-0.01 


0.10 


6 


-0.13 


0.11 


14 


RX J1546.6-3618 


-0.04 


-0.05 


0.07 


9 


-0.11 


0.10 


18 


RX J1547.6-4018 


-0.06 


-0.03 


0.09 


8 


-0.15 


0.12 


12 


Rho-Oph 
















RX J1620. 7-2348 


-0.08 


-0.04 


0.16 


6 


-0.12 


0.12 


17 


ONC 
















JW147 


-0.06 


-0.07 


0.12 


9 


-0.18 


0.07 


21 


JW589 


-0.12 


-0.20 


0.12 


4 


-0.18 


0.09 


11 


JW706 


-0.17 


-0.07 


0.13 


6 


-0.16 


0.08 


13 


Tau 
















1RXSJ052146.7 


-0.07 


0.01 


0.08 


7 


-0.12 


0.13 


11 


1RXSJ053650.0 


-0.18 


-0.14 


0.11 


4 


-0.23 


0.13 


11 


1RXSJ053931.0 


0.05 


0.05 


0.09 


10 


-0.03 


0.09 


13 



tra were obtained in the High-Efficiency mode (R~40 000) , 
positioning one of the fibres on the object and the other on 
the sky. Background light subtraction was performed us- 
ing the sky "spectrum" , since the instrument data reduc- 
tion software does not allow for the subtraction of inter- 
order scattered light. 

For all data sets, spectral extraction and data reduc- 
tion was achieved using the available instrument pipelines 
supplied by the observatories. Wavelength calibration was 
done using the spectrum of a thorium-argon lamp. The 
spectra were then Doppler corrected and set at a zero ve- 
locity. When more than one spectrum of the same star was 
available, they were co-added. Final spectra have S/N ra- 
tios between 100 and 200, and are of sufficient quality to 
derive stellar parameters and chemical abundances with a 
reasonable accuracy. In Tabled] we list fundamental phys- 
ical properties for all the observed stars in each target SFR. 

4. Spectroscopic analysis 

4.1. Stellar parameters and chemical abundances 



As in James et al 



stellar parameters and iron 
abundances for the targets were derived in LTE, using a 
grid o f plane-parallel, ATLAS9 model atmospheres (jKurucz 
19931) and the 2002 v ersion of the radiative transfer code 
MOOG ( SnedenI 1973). The meth odology used is described 
in detail in ISantos et al.l ( 20041 and references therein). 
The full spectroscopic analysis is based on the EWs of a 
set of Fe i and Fe ii weak lines, by imposing ionization and 
excitation equilibrium, as well as a zero slope between the 



abundances given by individual lines and their equivalent 
width. The errors in the stellar parameters we re derived 
using the same methodology as described in I Gonzalez 
(|1998[ ). 

For FGK dwarfs, the stellar parameters obtained using 
this methodology were shown to be compatible with other 
estimates in the literature. In particular, the effective tem- 
peratures derived are very close to the ones obtained us- 
ing recent applications of the infra-red flux method (e.g. 



Casagrande et al.ll2006l ) 



The results of these analysis are presented in Tabled] 
where we also compare the values derived using uvES spec- 

(|2006h sam- 



tra for stars in common with the James et al 



pie. We note that in their paper, James et al. used lower 
quality spectra taken with feros. A comparison of their 
results and our new ones shows that besides the larger 
error bars, the two datasets are in close agreement (well 
within error bars). For the duration of this paper, we will 
only consider the higher quality values derived from uvES 
spectra, for those stars having both feros and UVES mea- 
surements. 

A s was discovere d in the studies bv ljames et al. ( 2006 ) 
and (|Padgettlll99"i) . the derived micro-turbulence veloc- 
ities are considerably above the ones found for main- 
sequence dwarfs with similar temperature detailed in the 
literature. The cause of these discrepancies is not yet clear 
to us, although the high values may be reflecting the effect 
of magnetic fields ( Steenbock fc Holwegeij[l98ll i. 

Abundances for silicon (alpha-element) and nickel 
(iron-peak element) were derived using the method ology 
and atomic line lists described in iGilli et al.l (|2006h and 
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Table 3. Weighted average and dispersion of the metalhcities and [X/Fe] (X=Si, Ni) abundance ratios derived for each 
SFR. N and N' are the number of stars used for the determination of <[Fe/H]> and <[X/H]>, respectively. For Rho- 
Oph, since only one star was measured, the RMS of the values correspond to the rms of the abundance determination 
for the star. 



SFR <[Fg/H]> rms N <[Si/Fe]> a([Si/Fe]) <[Ni/Fe]> a([Ni/Fe]) N' 

cha -0.11 an 4 003 aoi -0.05 002 T 

CrA -0.06 0.05 3 0.01 0.01 -0.09 0.04 2 

Lup -0.05 0.01 5 0.02 0.02 -0.10 0.02 4 

Rho-Oph -0.08 0.12 1 0.04 0.16 -0.04 0.12 1 

ONC -0.13 0.06 3 0.00 0.09 -0.06 0.07 3 

Tau -0.05 0.11 3 0.04 0.04 -0.07 0.02 3 



Santos et al. (|2006l) . and are listed in Table[2j Analysis 
of these two elements was performed because they have 
several well defined lines in the spectra of solar type 
stars, and because they do not present any strong de- 
pendei ice of derived abundance with effectiv e temper- 
ature (|Bodaghee et al.l 120031: iGiUi et al.|[2006f) . This in- 
creases our confidence of uniformity in the comparison 
with abu ndances derived for solar- type field stars, pre- 
( 200^ . 



Table 4. Veiling determinations for the stars listed in 
Table[T]that were observed with UVES, flames or SOPHIE 
spectrographs. 



4.2. Veiling 

The spectroscopic analysis described above relies on the 
measurement of EWs for a set of Fe i and Fe il lines. 
This means that any contribution of veiling, related to 
stellar accretion phenomena, could lead to erroneous val- 
ues for the stellar parameters and effective temperatures. 
Although high levels of veiling are not ex pected for the 
WTTS in our sample ( Hartigan et al. 1995[ ). we err on the 
side of caution, and have performed a more detailed anal- 
ysis. 

The veiling (defined as the ratio between the excess 
continuum flux and the photospheric continuum contribu- 
tion to the spectrum) was derived usi n g a rn ethod based 
on that described in iHartigan et al. ( 1989h . With this 
method, the line depth of photospheric absorption lines 
are compared with a reference spectrum within a small 
wavelength band (a few tens of A wide) , and it is assumed 
that the veiling be constant in this region. The wavelength 
coverage of our spectra is almost complete from 4800 to 
7000A with a smah gap between 5700 and 5850A. For each 
object, a reference star (with similar temperature as the 
object) was rotationally broadened to match the v sin i of 
the object, and the veiling computed in intervals of 20A. 
The reference stars are dwarfs of similar effective tem- 
perature and metallicity to the targets, and were the ones 
used in previous chemical abund ances studies of stars with 
and without extra -solar planets (jSantos et al.ll2004 , 120051: 
Sousa et al.1l2006[ ). 



No clear trend in the veiling variation with wavelength 
was found for any of the stars. For all the objects, the 
veiling is centred near zero with some scatter (Fig. [2]). 
The values displayed in Table|4] represent the mean veil- 
ing bluewards of 5700A and redwards of 5850A. The error 
bars represent the standard deviation in each region. The 



Star 


Veil 


Veil 




<5700A 


>5850A 


Cha 






RX J1140.3-8321 


-0.01±0.04 


-0.07±0.06 


RX J1158.5-7754A 


0.03±0.04 


0.00±0.05 


RX J1159.7-7601 


-0.01±0.03 


-0.05±0.06 


RX J1233.5-7523 


0.02±0.01 


0.04±0.02 


CrA 






CrAPMS-4SE 


0.00±0.04 


-0.04±0.04 


RX J1839.0-3726 


0.00±0.05 


-0.05±0.05 


Lup 






RX J1507.2-3505 


0.00±0.05 


-0.03±0.05 


RX J1546.6-3618 


-0.01±0.04 


-0.05±0.05 


RX J1547.6-4018 


0.00±0.03 


-0.01±0.04 


RX J1523.4-4055 


-0.02±0.06 


-0.07±0.07 


Rho-Ohp 






RX J1620. 7-2348 


0.01±0.04 


-0.04±0.06 


ONC 






JW157 


0.00±0.05 


-0.06±0.08 


JW589 


O.OliO.lO 


-0.04±0.08 


JW706 


0.02±0.03 


-0.03±0.06 


Tau 






1RXSJ052146.7 


0.00±0.04 


-0.03±0.06 


1RXSJ053650.0 


0.00±0.03 


-0.02±0.05 


1RXSJ053931.0 


0.04±0.05 


0.00±0.05 



negative values and the highest error bars found in the 
red region are due to shallower photospheric lines and the 
presence of telluric lines in the object spectra, which are 
poorly fitted by the spectra of the rotationally broadened 
reference stars. 



5. Results and discussion 

5.1. Systematic errors or an abundance dispersion 
inside the sfr 

As can be seen from TablesH] and [H the general disper- 
sion observed in the [Fe/H] values of the stars in each SFR 
is quite small. The same is true for the abundances of 
Si and Ni (Tableland Fig. [3]). Possible exception is the 
Chamaeleon SFR, where the dispersion of the measured 
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5000 5500 6000 6500 

Wavelength (A) 

Fig. 2. Dependence of veiling with wavelength for five 
stars in our sample with different temperatures. The 
dashed and dotted lines show, respectively, the veiling at 
levels ±0.05 and ±0.10. Negative veiling occurs when the 
photospheric lines in the reference star are shallower than 
in the target. 



metallicities is quite high for the 4 stars studied (although 
still of the order of magnitude of the abundance uncer- 
tainties). Moreover, the abundances of the three species 
studied seem to be correlated, in the sense that higher iron 
abundances also correspond to higher [Si/H] and [Ni/H] 
values. 

A look at the temperature and metallicity values in 
Table[T] suggests, however, that a correlation may exist 
between metallicity and effective temperature. Using all 19 
stars in our sample, we find a slope of '^0.15dex/1000K 
for the points in a temperature-metallicity plot (in this 
fit, [Fe/H]=0.0 would occur for Tcs -5500 K). Although it 
may be unwise to use stars from different SFRs to conclude, 
this trend may suggest the existence of systematic errors 
in the derived stellar parameters. 

For a fixed metallicity, lower temperature stars tend 
to have deeper spectral lines. The strongest of these may 
have extended damping wings. The fitting of gaussian pro- 



files to measure the line-EWs could thus under-estimate 
the measurements, leading to lower metallicity estimates. 
Although this hypothesis could in principle explain the ob- 
served correlation, a look at individual lines in our spectra 
(Fig.[T]) shows that in the range of EWs measured, the ob- 
served spectral lines are well approximated by gaussian 
functions. This is actually expected, given that the large 
majority of the lines used have EWs below 100-150 mA. 

Furthermore, no correlation is seen in the upper enve- 
lope of the temperature- metallicity plo t of th e large sam- 
ple of stars presented in Santos et al.l (|2004h . whose pa- 
rameters were obtained using the same line-lists and model 
atmospheres. We can add that the metallicities and efec- 
tive temperatures derived for main sequence dwarfs using 
the methodology presented here are in good agreement 
with those found in the literature, d erived using other 
methods and model atm ospheres (e.g. ISantos et al.ll2004 : 



Casagrande et al.l l2006 



Morel et al.l (|2004) studied the detailed effects of stel- 



lar spots affecting the determination of stellar parame- 
ters and metallicity for RS CVn stars. Their results show 
that the final abundances can be affected, although not 
always strongly, by the presence of photospheric spots. 
This issue may have particular importance when study- 
ing young, magnetically active stars like the ones in our 
sample. These effects possibly contribute to the observed 
metallicity-temperature correlation. 

We have also verified the possibility that a small resid- 
ual veiling is causing the observed differences. For this 
purpose, we performed a simple test by increasing 3% the 
measured EW of iron lines used to derive stellar parame- 
ters for the two most metal deficient stars in the the Cha 
SFR (RX J1140.3-8321, RX J1158.5-7754a). The value of 
3% is actually a conservative upper limit for the observed 
veiling errors (note that the uncertainties listed in TableH] 
represent the observed rms in the different veiling points, 
and not the error on the average, which is much smaller). 
The resulting stellar parameters show only very small vari- 
ations (below 10 K in temperature and 0.01 dex in surface 
gravity), and their iron abundances increased by a mere 
0.02 dex. Clearly, small amounts of un-accounted veiling 
cannot explain the observed abundance dispersion. 

In any case, current data suggests that the observed 
dispersions may be the result of some un-identified sys- 
tematic effect in our data, and not to the existence of 
a real dispersion in the metallicities inside each SFR. This 
conclusion also implies that the observed dispersions must 
be seen as upper limits. 

5.2. Other elements 

Chemical abundances of several elements are useful 
in distinguis hing between different popu l ations in the 



Gala, xv (e.g. iNissenI 1200' 



2005 



Fulbright et al 



_'uhrmannll2004t iBensbv et al. 
Brewer fc CarnevI l200fil : iRich fc Origfial l2~ : 
20051) . erved differences are meant 



to reflect different star formation histories or stellar initial 
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-0.3 
[Fe/H] 



Fig. 3. Abundance ratios [Si/Fe] and [Ni/Fe] as a function 
of [Fe/H] for the six SFR (large circles) when compared 
with the same ratios for planet-host sta rs (filled triangles ) 
and "single" stars (open triangles) from lGilli et al. ( 20061 ). 
The dotted lines mark the solar abundances. 



mass functions in the different Galactic systems (|Carigi 
2006h . It is thus interesting to check whether the abun- 
dances of the elements studied in the current paper for 
the six SFRs are different from those found in field disk 
stars. 

In Fig.[3]we compare the abundance ratios [Si/Fe] and 
[Ni/Fe] listed in Table[3]with those found for f i eld pl anet- 
host stars and "single" stars from Gilli et al. I (|200d) . We 
remind the reader that Gilli et al. used the same proce- 
dure used in the current paper to derive the abundances 
of these elements (same model atmospheres, atomic line 
lists, etc), rendering such a comparison free from impor- 
tant systematic errors. 

The results show that no clear differences exist between 
the abundance ratios of the studied SFRs and those of field 
stars of the same metallicity. The chemical abundances in 
nearby SFRs seem to reflect those of typical, thin disk stars. 

5.3. The metallicity distribution of nearby SFKs 

In Table[31 we present the weighted average metallici- 
ties and abundance ratios for each SFR, derived from the 
values listed in Table[TJ In all cases the average abun- 
dances found are below solar. The metallicities for the 
six SFRS vary from -0.13 (ONC) to -0.05 (Lupus and 
Taurus). Intriguingly, this result suggests that metal-rich 
star-forming clouds in the solar neighbourhood are rare. 
The implications of this result are discussed in Sects. [?7?^ 
and l5.3.3l 



In Fig.|4]we compare the metallicity distribution of our 
sample of six SFRs (dashed line) with the same distribu- 
tion for 1250 stars i n the volurne-lim i ted CORALIE planet - 
search sample (see lUdrv et"aL l l2000t ISantos et al.l l200j ). 
The two distributions have similar average metallicities 
(—0.08 for the SFR sample, and —0.10 for the field stars). 
However, the metallicity distribution of the stars in the 
solar neighbourhood is clearly wider (cr=0.24dex), when 
compared with the one of the SFRs ((T=0.033dex]3. 

To have some insight into the statistical meaning of 
the observed difference, we performed two different tests. 
First, a simple Kolmogorov-Smirnov test was done to ver- 
ify the possibility that both samples belong to the same 
population. The result shows that P(KS)^14%, a value 
that is clearly not significant, likely due to the small num- 
ber of the observed SFRs. 

To better settle these issues we then performed a series 
of Monte-Carlo simulations. These were done to test the 
probability that in any given sample of randomly selected 
six metallicity points (taken from the field star sample) 
we would have a distribution similar to the one observed 
for the SFRS. In the simulation, the metallicity values were 
computed taking into account the observed metallicity dis- 
tribution for field stars presented in Fig.m T he method- 



ology used is similar to the one presented in lUdrv et al 
(|2003[ ). 



A total of 100 000 data sets were simulated, and for 
each of them we computed the values for CT([Fe/H]) of the 
six simulated points. We then compared the number of 
cases where the simulated (j([Fe/H]) is smaller than the 
observed one (0.033 dex). 

The results of these simulations show that the very 
small [Fe/H] dispersion observed is clearly at odds with 
the expected metallicity distribution of solar-type field 
stars. In our simulations, in only 0.02% of the cases we 
obtained a value of a([Fe/H]) lower than the observed 
0.033 dex (Fig.[5|). In other words, the very small disper- 
sion observed for the metallicity distribution of the ob- 
served SFRS is clearly significant. 

5.3.1. The small metallicity dispersion 

The small metallicity dispersion observed for the six star- 
forming re gions (0.033 dex) agrees w ith other recent stud- 
ies done bv lNieva fc Przvbillal ( 20061 ) , who derived similar 
values for the dispersion of the abundances in nearby OB 
stars (0.04 dex). When compared with local field dwarfs, 
a small dispersion is also observed in t he metallicities 
of open star clusters ( Twarog et al.|[l997 ). nearby young 
FGK stars and neb ular abundances (see compilation by 
Sofia fc Mev"Mll200lh . 

Together with the literature results, the data presented 
in this paper suggests that the material in the nearby 



^ The average metallicity and dispersion of the field star sam- 
ple is the same as the one found in the metallicity distribution 
for the subset of 94 CORA LIE sample stars wit h high precision 
metallicity estimates from lSantos et al] |200l). 
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[Fe/H] o{[Fe/H]) 



Fig. 4. Metallicity distribution of solar-neighbourhood 
stars in a volume-hmited sample (filled line - see text) and 
of the metallicities of the six SFRs analyzed in this paper 
(dashed line). In inset we show the cumulative functions 
of the two distributions. 



Fig. 5. Distribution of 100000 simulated [Fe/H] disper- 
sions (rms). The dotted line denotes the observed RMS 
of the 6 star-forming regions. In only 0.02% of the 
cases we obtained a dispersion below the observed value 
(0.033 dex). 



inter-stellar medium [ism] is very uniform. In turn, this 
suggests that at any given Galacto-centric radius, the iSM 
is likely to be well mixed, and that mixing timcscales are 
very short when compared with the time betwee n dijferent 
enrichment events (see discussion in IRov fc Kunthlll99,4 
Wielen et al. 19961 ). At a given epoch and Galacto-centric 
radius, the chemical abundances in the iSM are likely to 
be uniform. 

If we consider that nearby SFRs are good representa- 
tives of the current metallicity in the solar circle, the small 
scatter and slightly lower than solar abundances observed 
further imply that the metallicity distribution of main- 
sequence stars in the solar neighbourhood (with a much 
larger metalli city dispersion) is st rongly affected by stellar 
diffusion (e.g. Wielen et al. 1996h . 



5.3.2. The sub-solar metallicity of sfrs 

Although the average metallicity of the observed SFRs 
is below solar, in Sect. 15. II we have seen that some de- 
pendence seems to exist between the measured metallic- 
ities and effective temperatures (a positive gradient of 
~0. 15 dex/ 1000 K is observed). Since our stars are all 
cooler than the Sun (4664<Toff <5524K), this could im- 
ply that some under-estimate of the metallicities is done 
in our parameter estimate. We cannot exclude that the 
average metallicity of the observed SFRs may be closer to 
solar, although we can reasonably exclude a value that is 
strongly above. 



Our results 
support) some 



thus do not exclude (and somehow 
models of Galacto- chem ical evolution 
2005f ) that predict 



(jChiappini et al.l 120031 : iCarigi et al 



that a slow chemical enrichment has happened in the 
Galactic disk during the last 4.5 Gyr. We should add, how- 
ever, that models of chemical evolution of the galaxy often 
talk about oxygen as metallicity proxy, while in the cur- 
rent paper we measured iron abundances. Contrarily to 
oxygen, iron is formed in SN la events, that are supposed 
to continue even once the star-forming rate has slowed 
down. It will thus be interesting to explore the impact of 
our results on the models of galactic chemical evolution. 

If, on the other hand, the sub-solar abundances ob- 
served are real (or at least if the abundances are lower 
than those expected by models), the results could fit well 
a scenario where the Sun was formed in an inner region of 
the M ilky- Way's disk (later on experiencing orbi tal diffu- 
sion - I Wielen et al. 19961 : Edvardsson et al. 1993), or that 
the nearby iSM suffered from a recent in-fall of metal-poor 
gafl Sub-solar abundances wer e recently found for oxy- 
gen and carbon in B- type stars (iNieva fc PrzYbillal l2006h 



[zyb: 

and in the local bubble (jOliveira et al.l 120051 ). although 



in the latter case dust depletion may account for the dif- 
ference. These sub-solar metallicities are however at odds 
with so me other studies of chemi cal abundances in H ii 
regions (|Esteban et al.ll2004 [ioosh . 



Although this latter hypothesis seems unlikely given th e 
large quantities of printine gas needed fsee lCarigi et al.|[2006l ) 
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Recent values for the solar oxygen abundance (lower 
than previously considered) also suggest that the ism 
abu ndances may be closer to solar than p reviously thought 
(e.g. IChiappini eralll2003t lAsplund et al. 2004). A recent 
review on this question is addressed by lAllende Prietol 
([20071). However, if we consider that heavy el ement set- 
tling has occurred in the Sun (Lodders..2003.: Bahcall et al.l 
2QM) decreasing its abundance with respect to the initial 
value, the lower-abundances problem can be even more 
important than suggested by the mentioned values. 

5.3.3. On the origin of metal-rich stars 

As mentioned above, current models of galactic chemical 
evolution predict that a slow chemical enrichment has hap- 
pened in the galactic disk during the past 4.5 Gyr. Indeed, 
the models cannot easily explain the existence of super 
metal-rich stars in the solar neighbourhood. 

Our results thus support the idea that super metal- 
rich stars (with metallicity clearly above solar) found in 
the solar vicinity may have been f ormed in inner Galactic 
disk regions. In a recent paper, Ecuvillon et al.l (2007) 
present evidence that nearby metal-rich stars have kine- 
matics si milar to those of the Hyades (metal-rich) stream 
(see e.g. Famaev et al.l EoOtI ) . This is also the case of 
stars with planetary-mass companions, that are found 
to be metal-rich when compared with s in gle field stars 
Santos et alll2004l: lEcuvillon et al.ll2007h . iFamaev et al 



2007D and lEcuvillon et al.1 (|2007l ) interpret these results 



on the basis of a scenario in which metal-rich stars in the 
solar neighbourhood (and in particular planet-hosts) could 
have formed in a more metal rich ISM and then suffered 
radial displacements due to a non-axisymmetric compo- 
nent of the Galactic pote ntial, pushing them into the solar 
neighbourhood (see also I GrenonI 1 1 9991 ) . 

In this sense, it is interesting to see that for many 
chemical elements the abundance ratios [X/Fe] seem to 
present a rising trend as a function of [Fe/H] in the metal- 
rich domain, i. e., for [Fe/H]>0.0 (jBensbv et al.l l2005t 



Gilli et al. 2006h . Such trends seem very difficult to ex- 



plain with cu rrent models of galactic chemical evolution 
( Carigi 20061 ). Curiously, for a given [Fe/H], higher abun- 
dance ratios are als o observed for s t ars in other galac- 
tic populations (e.g. iFuhrmann 2004 IZoccali et al. I I2OO6I: 
Rich fc Origiiall200,^ . 

Also interestingly, there are suggestions that lithium 
abundances in planet-hosting, metal-rich stars show a sys- 
tematic difference with respect to field st ars in the tem- 



perature range between 5 600 and 5850 K (jlsraelian et al 



2004HChen fc Zhaolliooeh . We can speculate that the ex- 
planation for this difference can be related to a different 
lithium evolution in the inner Galactic disk. 

6. Concluding remarks 

In this paper, we present stellar metallicities for six 
nearby SFRs, derived from high resolution spectroscopy of 
young WTTSs. The results show that all the studied SFRs 



have metallicities between —0.13 and — 0.05dex, with a 
very small scatter of 0.033 dex around the mean value of 
<[Fe/H]>=:— O.OSdex (compatible to solar within the un- 
certainties) . 

The analysis of the abundances of silicon (alpha- 
element) and nickel (iron-peak element) also suggests that 
the abundance ratios in the studied SFRs are typical of 
those found in solar neighbourhood thin disk stars of sim- 
ilar metallicity. 

Together with other results from the literature, these 
observations suggest that the chemical abundances in the 
nearby ism are extremely uniform and not strongly above 
solar. Such a conclusion may have important implications 
for the models of Galactic evolution. 

The results presented here have also important impli- 
cations on planet searches around young (active) stars. 
Although these searches are se verely limited in cur 
rent radial-velocit" 



Santos et al. 200 



ity 
i 



(e.g. 'Saar fc Donahud 1X997 



surveys 

Paulson et al., ,2002 ' 3, several instru- 



ments are now being built that will change this situa- 
tion. In particular, new generation adaptive optics sys- 
tems like the ESO sphere planet finder or the Gemini 
GPI project, as well as a new generation of interferometric 
instruments like prima (eso), will allow us to study, in 
unprecedented detail, the existence of long period planets 
around young solar- type stars. The rarity of metal- rich 
nearby star-forming regions may limit the goals of these 
projects if the mctallicity-giant planet connection is still 
present for systems with long orbital periods. 

In this context it will also be very interesting to study 
in detail t he abundances of nearby young, post-T Taur i 
stars (e.g. IZuckerman fc Sond [200I iTorres et aLll2006h . 
Such studies would also help to set tle the question abou t 
the metallicity of the local ism fe.g. ISofia fc MevenbOOll ). 

If stars with planets have originated from inner galac- 
tic regions, it may be interesting to understand how plan- 
etary evolution (and survival) may suffer with stellar or- 
bital diffusion. This conclusion may further have impact 
on the understanding of the frequency of planets in the 
Galaxy. 

As a very important by-product, the observations pre- 
sented in this paper can further be used to investigate the 
effects of the metallicity on the properties of stellar popu- 
lations such as rotation, multiplicity degree, or magnetic 
activity. These studies will be presented in a separate pa- 
per, James et al. (in prep.). 
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See however recent detection by ISetiawan et al.1 (|2008l ) 
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